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Abstract 

In this paper, we use the fractional calculus to discuss the fractional mechanics, where the time 
derivative is replaced with the fractional derivative of order v. We deal with the motion of a body 
in a resisting medium where the retarding force is assumed to be proportional to the fractional 
velocity which is obtained by acting the fractional derivative on the position. The fractional 
harmonic oscillator problem, fractional damped oscillator problem and fractional forced oscillator 
problem are also studied. 



* Electronic address: |mimip4444@haninail.net| 
^Electronic address: Isosostl004@nate. coral 



1 



I. INTRODUCTION 



Fractional calculus implies the calculus of the differentiation and integration whose order 
is given by a fractional number. The history of the fractional derivatives goes back to the 
seventeenth century. There are good textbooks for the fractional calculus [1-7] . 

During about thirty years or so, fractional calculus has attracted much attention due to 
its application in various fields of science, engineering , finance and optimal problem. For 
various applications of fractional calculus in physics, see [8-18] and references therein. 

Many application of fractional calulus amount to replacing the time derivative in an 
evolution equation with a derivative of a fractional order. One of the problems encountered 
in the field is what kind of fractional derivative will replace the integer derivative for a 
given problem. Non-conservative Lagrangian and Hamiltonian mechanics were investigated 
by Riewe within fractional calculus [19,20]. Besides, Lagrangian and Hamiltonian fractional 
sequential mechanics, the models with symmetric fractional derivative were studied in [21,22] 
and the properties of fractional differential forms were introduced [23]. 

Two types of fractional derivatives , namely Riemann-Louville and Caputo, are famous. 
Mathematicians prefer Riemann-Louville fractional derivative because it is amenable to lots 
of mathematical manipulations. However, the Riemann-Louville fractional derivative of a 
constant is not zero, and it requires fractional initial conditions which are not generally 
specified. In contrast, Caputo derivative of a constant is zero, and a fractional differential 
equation expressed in terms of Caputo fractional derivative requires standard boundary 
condition. For this reason, physicists and engineers prefer Caputo fractional derivative. 

Recently, an extension of the simplest fractional variational problem and the fractional 
variational problem of Lagrange with constraints was obtained by using Caputo fractional 
derivative [10]. Even more recently, this approach is extended to Lagrangian formalism [24] 
with linear velocity. The fractional Hamiltonian formulations were presented for discrete and 
continuous systems whose dynamics were defined in terms of fractional derivative [25-31]. 

In this paper we will use the Caputo fractional derivative to discuss the fractional damped 
oscillator problem and the forced oscillator problem. The Caputo fractional derivative is 
defined as 




(1) 



where n = [i/] -|- 1 and [x] implies a Gauss symbol. 
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It can be easily checked that the fractional derivative satisfies the following : 

I (/3 = 0) 

Prom now on, we will restrict our concern to the case that < z/ < 1. Then the Caputo 
fractional derivative is given by 

II. INVERSE OPERATOR APPROACH 

The simplest first order differential equation is given by 

^ = /W (OO) (4) 

If we consider ^ — D], which means the fractional derivative with a — 1, the eq.(4) is 
re- written as 

DlN{t) = fit) (5) 
Prom the following relation for the fractional derivative 

D^D^^Dr^ D^,=Id, (6) 

we can solve the eq.(5): 

N{t)^D;'f{t)^ [fiOdC, (7) 

Jo 

where we insert a = — 1 in the eq.(l). 

Let us consider the first order fractional kinetic equation 

Nit)-NiO)^-c''DrN{t), (8) 

which is solved [32] as follows: 

N{t) = N{0)E,{-c''t'') (9) 
Here Ei,[x) is known as Mittag-Leffler function defined by 



k 

X 



^.W-Ef;;;;^ (10) 
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With a help of the fractional calculus, we are able to construct the classical mechanics 
related to the fractional calculus. Let us introduce the fractional velocity v{t)a,n.d fractional 
acceleration a{t) as follows : 

v{t)^D-^x{t), a{t)^Dy{t), (11) 

where x{t) is a fractional position and the time derivative is a Caputo's fractional deriva- 
tive. 

In the fractional mechanics, Newton's equation is defined by 

F^ma^ mD'^v, (12) 

where m is a mass of the body. When the force is constant, the body moves with the 
constant acceleration ^. It should be emphasized that the force is not necessarily constant 
and, indeed, it may consist of several distinct parts. Now let us consider the vertical motion 
of a body in a resisting medium in which there again exists a retarding force proportional 
to the fractional velocity. Let us consider that the body is projected downward with zero 
initial velocity v{0) = in a uniform gravitational field. The equation of motion is then 
given by 

mD^v — mg — kv (13) 
If we integrate the eq.(13), we obtain 

v{t) = gD;^{l) - ^D;^{v{t)) (14) 

lib 

If we multiply both sides of the eq.(14) by {—^YD^'^^ and sum up both sides for m from 
to oo, we obtain 

j2{--rDrMt) - E (--)'^+^A-^"'-'^(t) = ^ E(--)"A-^"'-'^'^(i), (15) 

n=0 ^ m=0 ^ n=0 ^ 



which yields 

m n r / i- M 

(16) 



1 - E,. -^f^ 



\ m J 

Here the asymptotic expansion of Ei,{—x) is given [33] by 



^ ^0 ^" 
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where bnii') is defined by 

l r-oo I CO 

— / dw[E2u{—'w^) COS kw + wE2u i+ui—w^) sin kw] — — 'S2bn{i-')k"' (18) 

with 

POO 

ho{^) = / E2,{-e)dt (19) 
Jo 



and 



The function defined by the eq.(20) is a generahzation of Mittag-Leffler function and first 
appeared in the work of Wimam [34] . 

The asymptotic expansion (17) shows us 

hm v{t) = (21) 

which is a terminal velocity in the fractional mechanics and gives a same result as that in 
the ordinary mechanics. 

Now let us introduce the fractional harmonic oscillator problem : 

mD^^xit) = -mw'^x{t), 0<u <1 (22) 

with the following initial condition 

x{0)^A, {D'^x){0)^vo (23) 

If we integrate the eq.(22), we obtain 

D'^x{t) -vo = -w'^D;''x{t) (24) 

If we integrate the eq.(24) again, wc get 

x(t) - x(0) = voD-^il) - w'^Df''x{t) (25) 

If we multiply both sides of the eq.(25) by {-w'^)"' D^^""" , we get 

{-w^rD-^"'''x{t) - {-w^){-wTD;'^"'^''^''x{t) 
= A(-«;2)"^D-2"^'^(1) + i;o(-t(;')"^A"^''"^'^"(l) (26) 
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Now summing up both sides of the eq.(26) for m from to oo, it gives 

oo oo 

TO=0 m=0 

oo oo 

= ^ E {-w^'Df-il) + ^0 E {-wTD^^'^^^'^il), (27) 



m=0 m=0 

which yields 



x{t) = AC^wt") + '^S.iwe) (28) 



and Mittag-Leffler cosine and sine function is defined as 

c(x) = y \ ^ , ^ , sjx) = y / ; f — -- (29) 

Mittag-Leffler cosine and sine function are also written as 

a(x) = \[E,{ix) + E,{-ix)], S,{x) = j^ESx) - E,{-ix)] (30) 



III. FRACTIONAL DAMPED OSCILLATOR 

The damped oscillator is the simplest classical model of motion with dissipation which 
corresponds to friction force proportional to velocity. Analogously, we can consider the 
fractional damped oscillator defined by 

mD'l''x{t) = -mw^x{t) - 2771^0"^ x{t), 0<p <1 (31) 

with the following initial condition 

xiO)^A, {D'^x){0)^vo (32) 

Now we will solve the eq.(31) by using the series approach. Let us assume that the 
solution of the eq.(31) is given by 



oo 



x{t) = y c^r^' (33) 

n=0 

Inserting the eq.(33) into the eq.(31), we have the following recurrence relation: 

an+2 + 27a„+i + w^ttn = 0, (34) 

where 

On = c„r(l + nv) (35) 



The recurrence relation (34) is re- written as 

a„+2 - pan+i = g(a„+i - pa„) , (36) 
where p, q are two solutions of the quadratic equation 

+ 2-ft + w^ = (37) 

Thus three possibilities exist. 

Case I : over damped oscillator In this case 7^ > and we put 



p — —7 + Y^7^ — w^, q — —7 — yj^"^ — w"^ (38) 
Case II : critical damped oscillator In this case 7^ = w'^ and we put 

P = g = -7 (39) 
Case III : under damped oscillator In this case 7^ < and we put 



p = —J + iy w"^ — g = — 7 — iyiw^ — 7^ (40) 
Iterating the recurrence relation yields 

an - p"ao = (ai - pao)p" ^ H - , (41) 



which gives the following solution 



r y4ip" - Biq'' {p ^ q) 

{ p'^ao + (ai - pao)np'^~^ {p = q) 



with 



(42) 



ai-pao ai - pao ..„s 

^1 = ao H , Bi = (43) 

p-q p-q 



Thus, the position operator has the following form : 

r Vo^E,{pt^) - --^E,{qt^) {p ^ q) 
x(t) = i ^--^ ^ ^ ^'^^^ (44) 

( AE^iptn + Bi^f^E^^^iptn {p = q) 

In order to discuss the under damped case, we introduce the following functions: 



00 00 



'I* > = „?. S r(i + (2„ + mMr(i + 2«)r(i + m)" * ' • '^^^ 
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S (a bin = V V ( l)"r(2 + 2n + m) m.2n+l.{2n+m+l)u /^gN 

Two functions satisfy the following relation: 

a(a, -blf") = a(a, 5,(a, -b\f) = -^,(a, (47) 

and 

a(a, 6|0) = 1, S,{a, b\0) = 0, [^^^'.(a, 6|r)]t=o = a, pr-^.la, 6|r)]t=o = b (48) 
Using these function, for the under damped case, the position is given by 

x{t) = AC4-^,w^\n + ^^^^S^-'y.wiin, (49) 

Wi 

where 

Wi — ^Jw^ — (50) 

When u goes to 1, we have 

a(a, 6|r) ^ e"* cos 6t, 5^(a, b\r) e"* sin 6t, (51) 
which implies that the eq.(49) is consistent with the classical result in this limit. 

IV. FORCED OSCILLATOR 

Now let us consider the forced oscillator with the driven force which is expressed in terms 
of Mittag-Leffler cosine function. For simplicity, let us restrict our discussion to the over 
damped oscillator case. Then the equation of motion is 

mD^^xit) = -mw^x{t) - 2m-fD'^x{t) + FoC^{(f)t''), 0<u<l (52) 

with the following initial condition 

x{0) = A, {D'^x){0)=vo (53) 

Let us assume that the solution of the eq.(52) is given by 



oo 



x(t) = 5:cX^ (54) 

n=0 



Inserting the eq.(54) into the eq.(52), we have the following recurrence relation: 



a2k+3 + 2702fc+2 + w^a2jt+i = 0, {k = 0, 1, 2, ■ ■ ■) 

«2it+2 + 27a2fe+i + w^a2jt = —(-0^)'', (A; = 0, 1, 2, • • •), 

m 

where 

a„ = c„r(l + nu) 
The recurrence relation (55) is re-written as 

0'2k+3 — PO'2k+2 — Q{(l2k+2 — P0'2k+l) , 

a2k+2-pa2k+i = q{a2k+i - pa2k) ^ (-0^)'', 

m 

where 7^ > and 



p — —7 + ^7^ — w^, g = —7 — ^J^^ — w^ 



Solving the recurrence relation (57) yields 



0'2k+i - P0'2k = ? (oi - pao) + 



m \ + ((P' 



a2fc+2 - pa2k+i ^ q (oi - pao) H . o , ^2 / 

my + 0^ y 

Inserting the first relation of eq.(59) into the second relation of eq.(59) yields 

a2.+2 = pW + {p + q) (a, - pa^ + - ^(/^ ^2) - 0^)(-0^) 

Solving the eq.(60) and inserting into the eq.(59), we have 

a2k = (ao + r/)/>2' + ^q^' + C(-0')' 
a2k+i = (aoP + r7')p'*^ + («l-P«o + O9'*^ + C'(-0')^ 



where 



1 ( qFo \ Fo ^ 2 



p — g \ m(g^ + 0^) / m(g^ + 0^)(p^ + 0^) 

C = (ai - pao + ] C = ^0 . _ ^2. 

q — p \ m{q'^ + 0^) / ' m(g2 + 0^) (p^ + 0^) 



m(g2 + 0^) ' m(g2 + 0^) 

Thus, the position has the following form : 

x{t) = {ao + v)CKipn + ^CK{qr) + CC,{<t>f) 
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+ (ao + + 



ai - pao + C 



(63) 



Q 



where ao and ai is determined from the initial condition and fractional hyperbolic cosine 
and sine function is defined as 



V. CONCLUSION 

In this paper, we used the fractional Caputo derivative to discuss the fractional mechanics, 
where the time derivative is replaced with the fractional derivative of order v. We dealt 
with the motion of a body in a resisting medium where the retarding force is assumed to be 
proportional to the fractional velocity which is obtained by acting the fractional derivative 
on the position. In this case the terminal velocity was shown to be the same as that in the 
ordinary mechanics. 

We used the inverse operator approach to solve the fractional harmonic oscillator problem. 
For over damped and critical damped case, the position was shown to be expressed in terms 
of Mittag-Leffler function. However, for the under damped case, it was shown that the 
position is expressed in terms of new functions named Ci,{a, bit"), Si,{a, bit") due to the fact 
that 



Finally, we used the series approach to solve the fractional forced oscillator problem with 
the driven force which is expressed in terms of Mittag-Leffler cosine function. 
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